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Abstract. Over the past three years, we have worked on developing a well-characterized 30 S 
radioactive beam to be used in a future experiment aiming to directly measure the 30 S(a,p) stellar 
I ' reaction rate within the Gamow window of Type I X-ray bursts. The importance of the S(a,p) 

reaction to X-ray bursts is discussed. Given the astrophysical motivation, the successful results 
^ ■ of and challenges involved in the production of a low-energy 30 S beam are detailed. Finally, an 

overview of our future plans regarding this on-going project are presented. 
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The 30 S(a,p) reaction is a significant link in the ap-process, which competes with the rp- 
process in Type I X-ray bursts (XRBs) [1], but the reaction rate is virtually unconstrained 
by experimental data. Not only are there no data directly measuring the 30 S(a,p) cross 
section in the literature (at any energy), but there are no conclusive experimental reports 
on the nuclear structure of the compound nucleus 34 Ar above the a-threshold with regard 
to any possible oc -resonances. Hydrodynamic models of XRBs indicate that variation of 
the theoretical reaction rate has significant consequences. 

XRBs are understood to result from thermonuclear runaway in the hydrogen- and 
helium-rich accreted envelopes on the electron-degenerate surfaces of neutron star bi- 
nary systems [2]. Accretion ensures a steady flow of fresh material will spread around 
the surface of the neutron star only to be buried by the continuous pile-up of more mat- 
ter. Equilibrium between a fierce gravitational pressure and the constant energy release 
of the /3 -limited CNO cycle is broken by a thin-shell instability, triggering the onset 
of explosive nucleosynthesis [3, 4] . Although powerful, these bursts do not disrupt the 
binary star system, hence X-ray bursters exhibit recurring episodes with hourly, daily, 
or more extended regularity, making them the "most common thermonuclear explosions 
in the universe" [5]. Computer simulations reproduce the energy release, burst profiles 
(rise time, peak wavelength, and decay curve), and the recurrence time-scales of these 
astrophysical phenomena. 



The theoretical 3 S(a,p) cross section at astrophysical energies is typically calculated 
using the statistical model of Hauser and Feshbach [6]. However, statistical models 
must assume an energy-dependent level density in the compound nucleus of interest, 
which is poorly constrained experimentally for the case of 34 Ar. In the case where a 
cross section is dominated by one or more narrow-resonant contribution(s), one ought to 
treat the problem via resonant reaction formalism and disfavor a statistical approach. 
Previous work indicates that for a-induced reactions on T z = ±1 (T z = (N — Z)/2) 
nuclei with A = 18, the cross sections are shown to be dominated by natural-parity, 
a-cluster resonances [7, 8, 9]. These a-capture trends are also observed at higher mass 
(22 < A < 30) on T z = 1 nuclei [10, 11], calling into question the accuracy of applying 
a statistical model to reaction rates such as 30 S(a,p). 

With nuclear reaction networks modeling the rp-process now extending beyond A = 
100 and involving a myriad of nuclear transmutations, one may wonder the extent 
to which variation of individual nuclear reaction rates within their uncertainties may 
give rise to gross systematic effects. By varying only the 30 S(a,p) reaction rate by a 
factor of 100, one model shows a different XRB profile, possibly accounting for the 
double-peaked structure in the bolometric luminosity of some rare systems [12]. The 
30 S(a,p) reaction is also found to alter the crustal composition of neutron stars [13], 
which influences the compositional inertia of recurrent XRBs [14]. A parametric study 
on the nuclear inputs to XRBs quantifies the effects of varying many reaction rates, 
indicating that 30 S(a,p) is one of the nine most influential nuclear processes occurring 
below 56 Ni and one of eight processes below 56 Ni affecting the energy output by more 
than 5% [15]. Clearly, experimental measurements to constrain the 30 S(a,p) reaction 
rate are warranted. 



3U S BEAM PRODUCTION 

The physical chemistry of sulfur makes it extremely difficult to extract quickly and 
efficiently for re- acceleration as a secondary beam in ISOL facilities, and a method for 
producing low-energy (< 4 MeV/u) secondary beams well-characterized in phase- space 
has yet to be realized using high-energy fragmentation techniques. As detailed presently, 
the low-energy Center for Nuclear Study (CNS) radioactive ion beam (CRIB) separator 
facility of the University of Tokyo [16, 17] and located at the Nishina Center of RIKEN 
is capable of producing a 30 S RI beam suitable for studying the astrophysical 30 S(a,p) 
reaction. 

Heavy ion beams of a reasonably high intensity are extracted from the CNS 14 GHz 
HyperECR ion source [18], accelerated to an energy < 12 MeV/u by the AVF cyclotron 
(K = 70) at the RIKEN Accelerator Research Facility, and delivered to the entrance 
focal plane of CRIB (denoted 'F0,' see below). For the present work, we bombard a 
cryogenically-cooled 3 He gas target with a stable 28 Si beam, and via the 3 He( 28 Si, 30 S,)n 
reaction, we produce the RI species of interest. The target 3 He gas is confined by 2.5 /im 
Havar windows and force-flowed through a LN2 cooling system to an effective target 
temperature around 80-90 K [19]; the resulting 3 He gas density is a factor of ~3 greater 
than at ambient laboratory temperatures and stable against any density reduction effects 
induced by the energy deposition of the primary beam. 




The cocktail beam emerging from the production target is mainly characterized and 
purified in the experiment hall by two magnetic dipoles and a Wien (velocity) filter, with 
beam-focusing magnetic multipoles surrounding these elements (see Figure 1). These 
beam-line components are separated by four focal planes of interest. The primary beam 
focal point and the production target are located at FO, the dispersive focal plane between 
the two magnetic dipoles is denoted 'Fl,' the achromatic focal point after the second 
dipole 'F2,' and the location of the experiment scattering chamber after the Wien filter 
'F3.' Slits limiting the emittance and detectors for diagnosing the beam (such as silicon 
detectors and parallel plate avalanche counters (PPACs)) can be inserted or removed at 
all four and the latter three of the focal planes, respectively. The slits located at Fl are 
of particular interest as they clearly define the magnetic rigidity (Bp) of transmitted ions 
as well as their maximum momentum dispersion (Ap/p). 

As we conducted 30 S RI beam development tests in December 2006, May 2008, and 
July 2009 (two days each) varying many parameters to optimize results for 30 S, we 
will limit the discussion to highlights of the most noteworthy points. We tested three 
primary beams: 28 Si 9+ of 6.9 MeV/u at 100 pnA, 28 Si 10 + of 7.54 MeV/u at 10 pnA, 
and 28 Si 9+ of 7.4 MeV/u at 144 pnA, listed in chronological test order. All intensities 
quoted here were the maximum available at the time of the various tests, and due to 
beam-monitor counting limits at some focal planes, the actual intensity used may be 
lower where results are then normalized; previous work at CRIB indicates that, up to 
the rates of energy deposited in the production target for this work, linear beam-current 
normalization is applicable. 
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FIGURE 2. The plot shows the particle flight time in nanoseconds on the abscissa and the residual 
energy in million electron volts on the ordinate for various nuclear species in the cocktail beam at the 
achromatic focal plane F2. Although 30 S 16+ is clearly separated, the loci of other charge-states of 30 S are 
heavily contaminated. The dispersive momentum slits are set such that Ap/ p < 0.625%. 



We found that 3U S beam intensity shows a positive correlation with primary beam 
energy within this range, justifying our choice of the highest 28 Si beam energy available 
from the cyclotron for each test; the results of recent improvements bringing the practical 
AVF cyclotron K-value closer to its design value are evident comparing the two 28 Si 9+ 

beam energies listed above. Although dG jp for 3 He( 28 Si, 30 S) is negative within this 
energy range [20, 21], the combined energy dependence of the momentum straggling 
within the production target and the emerging charge-state distribution of 30 S more 
than compensate for the decreasing o(E); RI beam production and ion-optical transport 
simulations are consistent with this interpretation of our experimental data. 

The production target thickness was also optimized for 30 S yield. We experimentally 
found this value to be ~ 1.7 mg/cm 2 of 3 He, corresponding to a cryogenic gas pressure 
of 400 Torr. We tested pressures of 200, 300 and 400 Torr for 30 S production, and 
simulations as well as experimental tests with 3 He( 16 0, 18 Ne) indicate that significantly 
higher target thicknesses decreased RIB yield on-target. These results can be understood 
in the same framework as the aforementioned primary beam energy effects: despite the 
linearly increasing isotropic yield of 30 S for thicker targets, at some point the increasing 
momentum straggling and unfavorable shifting of the charge-state distribution from 
energy-loss considerations are dominate contributions to the overall intensity. 

Charge- state effects may seem unimportant since in principle we are free to tune 
the CRIB optics to any species of interest, and we may just choose the predominately 
populated ion-species of 30 S. To appreciate the importance of charge-state effects, one 
must consider that magnetic dipoles separate based on p/q, which, with a narrow 
momentum dispersion, effectively reduces to the charge-to-mass ratio A/q. Owing to 
an extremely high intensity and statistical straggling effects, various charge states of 
the primary beam can be found at nearly any momenta. Only the fully-stripped ion 
30 S 16+ (A/q ~ 1.875) is clearly separated from the leaky primary beam (see Figure 2), 



which can never have A/q < 2 ( 28 Si has N = Z and hence A = 2Z and q max = Z) 1 . 30 P 
arising from the 3 He( 28 Si, 30 P)p reaction, with a cross section of the same order as the 
3 He( 28 Si, 30 S) reaction [20], also shows up as an impurity for all but the fully- stripped 
case of 30 S. In December 2006, we determined that without measuring the energy or a 
significant energy-loss of the cocktail beam, it was impossible to fully separate 30 S 15 + 
from 28 Si 14+ at a satisfactory level to avoid false-positives. In May 2008 we could not 
purify 30 S+ 14 above the ~ 1% level even with use of the Wien filter. Although analysis of 
the July 2009 beam test is still underway, preliminary analysis indicates no significant 
improvement in the status of the 30 S 14+ beam purity (< 2%). Our results for 30 S 16 + 
have continued to improve each year, and we have successfully achieved ~ 10 4 particle 
Hz on target with ~ 25 — 30% purity with E beam = 30 ± 3 MeV. Unfortunately at these 
energies, the 16+ species of 30 S is very weakly populated (< 1%), accounting for its 
low intensity. 

As the production gas target has Havar windows, we considered ways of increasing the 
intensity of fully-stripped charge states of ions emerging from Havar foil. We measured 
the charge-state distribution of 28 Si beam ions in a thick carbon foil (550 /ig/cm 2 ) 
compared to Havar foil (2.2 /im) (Table 1); the energy-loss in Havar is slightly greater 
than carbon in this case, but the difference does not significantly effect our interpretation. 
For a 28 Si beam of 3.4 MeV/u (~ E beam of 30 S), it was found that transmission of highly 
charged states of 28 Si is improved through carbon foil compared to Havar foil with a ratio 
consistent with predictions of LISE++ (specifically, the Global+Leon model) [22, 23]. 
In our July 2009 test, we used a 2.5 /im Be foil after the production target, which when 
normalized for comparison with the May 2008 results, indicates an increase in the 30 S 16+ 
intensity by a factor of 2. Although one theoretically expects this intensity increase to 
be on the order of a factor of 10-20, the Be foil was partially broken and perhaps did 
not cover the full solid-angle of beam emission from the production target, possibly 
accounting for this deficiency. 



TABLE 1. Intensity of selected 
charge states of 28 Si after passing 
through Havar foil or carbon foil. 



Target 


Species 


Normalized 






pps @ 10 enA 


Havar 


28 Si 12+ 


1.075 x 10 8 


Havar 


28 Si 13+ 


6.013 x 10 7 


Havar 


28 S jl4+ 


3.901 x 10 6 


Carbon 


28 Si 12+ 


1.758 x 10 8 


Carbon 


28 Si 13+ 


1.300 x 10 8 


Carbon 


28 si 14+ 


4.365 x 10 7 



For simplicity, we quote A /q without units. 



PLANNED EXPERIMENTAL SETUP 



We plan to measure the 4 He( 30 S,p) reaction 2 concurrently with 30 S+a resonant elastic 
scattering using the thick target method [24]. The measurements will scan from the top 
of the 2 GK astrophysical Gamow window down to 1 GK, requiring a 30 S Et, eam = 32.3 
MeV on-target; while models indicate the explosive conditions in XRBs range from 
0.4 < T < 1.3 GK [25, 5], the predicted cross sections at the lowest energies are much 
too low to be measured with the presently developed 30 S beam. Indeed, Hauser-Feshbach 
calculations indicate that we require a 30 S beam with intensity 10 5 particle Hz to get 
reasonable (oc,p) reaction statistics in the 2 GK region over the course of an 11 day 
experiment. But, with the presently-developed 30 S beam at 10 4 particle Hz, we may 
reasonably measure the resonant a scattering of 30 S and at least put a meaningful upper 
limit on the (a,p) cross section; as the Fisker, Thielemann & Wiescher model shows 
significant effects on the overall XRB light-curve for 30 S(a,p) cross sections 10 2 above 
the Hauser-Feshbach rate [12], we may at least confirm or rule-out such a rate with the 
planned experiment. 

We will measure the reaction ejecta on an event-by-event basis in coincidence with 
30 S ions detected in two upstream beam monitors (PPACs or Multi-Channel Plates); as 
these beam monitors are limited to a counting rate of ~10 6 particle Hz [26], this puts an 
upper limit on the maximum beam intensity on-target. The original experiment proposal 
called for a semi-cylindrical 4 He gas cell and AE-E silicon telescopes. However, ejecta 
data collection in this setup is limited to a laboratory solid angle of ~0.15 sr, and 
owing to the poor timing resolution of silicon detectors, clearly identifying (a,po) 
events or distinguishing inelastic a events may be challenging. Our group is currently 
developing an active-target method for helium gas, adding Gas Electron Multiplier 
(GEM) capabilities to a newly designed Multi-Sampling and Tracking Proportional 
Counter (MSTPC) based on previous research [27]. The active-target fill-gas is 90% 4 He 
and 10% C0 2 . The GEM-MSTPC will allow for the full track-reconstruction of beam 
particles, reaction points, recoil nuclei, and ejecta, as well as full energy measurements 
of the ejecta in silicon detectors, allowing us to clearly identify reactions not transiting 
directly to ground-state levels. The silicon detectors in the GEM-MSTPC cover ~15% 
of the laboratory solid angle, greatly increasing our data collection statistics compared 
to the originally proposed setup. 

SUMMARY 

We successfully developed a 30 S RI beam of 10 4 particle Hz of ~25% purity and 
Ebeam = 30 ± 3 MeV. In the spring of 2010, we will conduct an experiment on the 
first measurement of the 30 S+ 4 He system, the results of which will be applied to the 
astrophysical 30 S(a,p) reaction rate. 



Although the experiment will be performed in inverse-kinematics, for nomenclatural convenience we 
may also write these reactions in normal kinematics (e.g. 30 S(a,p)). 
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